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ABSTRACT 


The objective of the research under this Grant has been the development 
of approaches which are alternatives to current design methods which rely 
heavily on linear quadratic and Riccati equation methods. The main alterna- 
tives were classified into two broad categories: 

Local Multivariable Frequency Domain Methods 
Global Nonlinear Optimal Methods 

The progress made toward the original research objectives within the 
one year grant period is described briefly in this report. More detailed 
information is available in the following works, which were essentially 
completed during the grant period. 

(1) T. C. Brennan and R. J. Leake, "Simplified Simulation Models for Con- 
trol Studies of Turbojet Engines," Notre Dame Electrical Engineering 
Technical Report No. EE- 757, 88 pages, November 1975. 

(2) J. C. Shearer, "An IBM 370/158 Installation and User's Guide to the 
DYNGEN Jet Engine Simulator," Notre Dame Electrical Engineering M.S.E.E. 
Thesis, 44 pages, November 1975. 

(3) R. J. Leake and J. G. Allen, "Optimal Regulators and Follow-Up Systems 
Determined from Input-Output Signal Monitoring," Notre Dame Electrical 
Engineering Memorandum No. 7403, 16 pages, October 1975. 

(4) R. R. Gejjl and M. K. Sain, "Polynomial Techniques Applied to Multi- 
variable Control of Jet Engines," Notre Dame Electrical Engineering 
Technical Report No. EE-761, 118 pages, March 1976. 

(5) V. Seshadri and M. K. Sain, "Multivariable System Compensation Includ- 
ing Characteristic Methods for Jet Engine Control," Notre Dame Electrical 
Engineering Technical Report No. EE-763, 114 pages. May 1976. 

(6) M. K. Sain, et. al., "Alternative Methods for the Design of Jet Engine 
Control Systems," Proceedings of the 1976 Joint Automatic Control Con- 
ference, Purdue University, pp. 133-142, July 1976. 
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ALTERNATIVES FOR JET ENGINE CONTROL 


a. Polynomial Techniques 

R. R. Gejji and M. K, Sain 

b. Characteristic Loci Techniques 
V« Seshadrl and M. K. Sain 


la. POLYNOMIAL TECHNIQUES 


Completion of the first year's work on Grant NSG-30A8, occurring 
February 29, 1976, was simultaneous with the completion of preliminary 
Investigations on polynomial techniques applied to jet engine control. 

This work was carried out by R. R. Gejji under the direction of M, K. Sain 
and served also as Mr. Gejji’ s Master of Science thesis activity. As a 

research activity, the effort was jointly supported by this grant and by 

« 

the National Science Foundation under Grant GK-37285. Effort under this 
grant was focused upon the development of a realistic, practical applica- 
tion for the polynomial design method in the area of jet engine control, 
while the NSF support was used to provide some flexibility for theoretical 
and algorithmic extensions of polynomial concepts as they became nc ,ded. 

Before presenting an abstract of this wo’k, it may be useful to provide 
a brief notion of the types of problems and methods involved. 

Basically, polynomial techniques have to do with the solution of trans- 
fer function equations of the form 

A(s)X(s) » B(s) 

where A(s) and B(s) are known matrices each of whose elements is a ratio of 
polynomials in s, with real coefficients. X(s) is an unknown matrix in 
the same class and is to be determined. For feedback control of jet engines, 
A(s) and B(s) arc determined from engine data around the operating point of 
interest and from specifications on the performance of the feedback control 
system; and X(s) embodies information about the control compensations required 
to satisfy the specifications. If it were merely desired to solve exhaustively 
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the linear equation above, it would be possible to take advantage of the field 
structure of ratios of polynomials and proceed according to the classical Intro- 
ductory methods of linear equations. Unfortunately, this Is not the case, since 
in control problems there is also V^y interest in 

(i) whether X(s) is physically realizable, 

(li) where the poles of X(s) are located in the complex plane, 

(ill) and, if more than one X(s) is possible, how many dynamical 
elements are required to construct X(s). 

These issues immediately remove the solution from an elementary classification. 
One approach which has developed in recent years to handle issues such as (i) , 
(ii), and (ill) is what we have been calling the polynomial technique. The 
idea can be simply illustrated by writing 

A(s) >= N^(s)/a(s), 

B(s) ** Ng(s)/b(s), 

X(s) = N(s)d"^(s) 

where N. , N , N and D are matrices of polynomials and a and b are polyno- 
A D 

mials. The representation for X(s) is called a matrix fraction and is quite 
helpful in answering advanced questions about X. With these conventions, the 
design equation becomes 

lb(s)N (s) I -a(s)N (s)]I N(s)l “ lO] 

* • ® L >>(=)J 

so that pairs of corresponding columns from N and D can be understood 
algebraically as solutions to a homogeneous equation having a polynomial coef- 
ficient matrix. 

The work discussed in the following abstract concerns the specification 
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of A(s) and B(s)i the calculation of homogeneous solutions to this polynomial 
equation, and the control theoretic Interpretation of these solutions. 

Abstract 

R. R. Gejjl and M. K. Sain, "Polynomial Techniques Applied to Multivariable 
Control of Jet Engines, Technical Report No. EE-761, Department of Electrical 
Engineering, University of Notre Dame, Notre Dame, Indiana 46556, March 22, 1976. 

One way to approach the design of linear multivariable control systems is 
by means of exact model matching, where desired closed loop performance Is 
expressed in terms of a specified transfer function matrix. One questions which 
is often raised about such an approach involved the practicality of specifica- 
tion of such a desired, closed loop performance matrix. Another, related, 
question concerns the possibility of determining the existence of compensators 
to achieve such a performance, as well as giving a finite cnumoratln of them 
and selecting those which arc, for example, minimal or stable. 

This study addresses itself to these points, in reverse order. Research 
by Sain in 1975 has established that determination of the existence of compen- 
sators and a finite enumeration of them is possible within the context of free 
modules over polynomial rings. Algorithms to accomplish this in theory are 
also available. The first purpose of this work was to construct straight- 
forward computer software to check the workability of these theoretical algo- 
rithms. This has been accomplished in both FORTRAN and PLl, with listings 
for the former included herein. Tt was found that these programs were success- 
ful on the sort of small problems which often appear in the literature. The 
second purpose of the work was to evaluate these first-generation software 

efforts on a realistic practical problem. This has been accomplished in the 
• • 

context of multivariable control of models derived from the FIDO turbofan 
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engine by Michael and Farrar. It was established that specification of 

desired closed loop performance matrices is a reasonable idea for this appll-* 

cation, and, in fact, that in at least c.ie case it has already been done. 

This means that exact model matching is not a vacuous Idea, at least for 

jet engines. Efforts were then upended to determine if the first generation 

software would compute compensators, which had been proven to exist. The 

results here were promising but inconclusive. Though the software would 

yield ans'Atrs of the degree crnected (a definitely nontrivial achievement), 

eigenvalue tests and other checks initiated to verify accuracy were decidedly 

%• 

disappointing. 

Two conclusions are possible: (1) that some oversight remains undetected 

in the software or (2) that straightforward computer software is unequal to 
this application. The second of these pcssibilities appears to be the more 
probable at this time, though the first cannot be eliminated with certainty. 
Studies are under way to develop more sophisticated software. 

In any event, the results in terms of problem formulation are promising 
enough, and the computer results close enough to suggest that further work 
on this problem would be worthwhile. 

Remark 

Contained in this report are two results which are ready for presentation 
to the technical community. The first result, on discovering and fully formu- 
lating what seems to be an important exact model matching problem in jet 
engine control, is being prepared for the Nineteenth Midwest Symposium on 
Circuits and Systems, which will take place in the summer of 1976. The second 
result, on development of FORTRAN and PLl software for this class of problems, 
is ocing prepared for the Fourteenth Allcrton Conference on Circuit and System 
Theory In the fall of 1970. 
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Remark 

The sorts of academic polynomial equations previously solved in the 
literature have, in the experience of the authors, taken at least an hour 
or two to solve by hand. Moreover, hand solution is tedious and error-proven. 
The fact that this software works successfully on such problems is, thus, 
nontrivial and represents an improvement of roughly two orders of magnitude 
in solution time for such cases. The necessity for additional software 
development for the jet engine example can be regarded as a temporary delay 
rather than as an essential obstacle. Further results in this regard are 
expected to be achieved over the summer. 

Remark 

There is also one result on polynomial methods in the next section. 


Ib. CHARACTERISTIC LOCI TECHNIQUES 


In a nanner alnllar to che studies on polynomial techniques t completion 
of the first year's grant work occurred almost simultaneously with the comple- 
tion of preliminary investigations on characteristic loci techniques applied 
to jet engine control. This work was carried out by V. Seshadrl under the 
direction of M. K. Sain and also served as Mr. Seshadrl* s Master of Science 
thesis activity. As a research activity, the effort was jointly supported 
by this grant and by the National Science Foundation under Grant GK-37285. 

Effort under this grant was focused on the application of the characteristic 
loci technique to realistic numerical data associated with jet engine models, 
while the NSF support allowed some basic Investigation into the nature of 
characteristic loci themselves. 

A word or two on the meaning of characteristic loci may be helpful before 
a reading of the abstract ^Ich follows. 

Recall the Nyqulst method for analysis of a plant having one input and 
ont output. The strategy is to study the behavior of the plant transfer func- 
tion g(s) as s travels an appropriately chosen Nyqulst contour roughly enclos- 
ing the right half of the s-plane. To gain appreciation for MacFarlane's 
generalization of Nyqulst’s Idea. It is convenient to observe that the trans- 
fer function g(s) actually defines the eigenvalue of the plant operator, as 
can be seen from 

[g(«)][l3 ■ g(8)(l], 

where [1] is the eigenvector of the operator [g(s) and G(s) is its eigenvalue. 

In generalizing to the case of a multivariable plant G(s). having an equal 
number of Inputs and outpus. MacFarlane determined to make Nyqulst plots of 
Che eigenvalues of G(s). These plots are called characteristic loci. 
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Abstract 

V. Seshadrl and M. K. Sain, '*Multivarlable System Compensation Including 
Characteristic Methods for Jet Engine Control, "Technical Report No. EE-763, 
Department of Electrical Engineering, University of Notre Dame, Notre Dame, 
Indiana 465^6, May 6, 1976. 

This study addresses two possible approaches to the problem of linear 
awltlvarlable control system design. The first approach belongs to the class 
of exact model matching designs; the second approach makes use of the charac- 
teristic methodologies of MacFarlane. 

In the exact mbdel matching case, we consider the problem of finding 
stable solutions G(s) to the transfer function luatrix equation 

Gj^(v)G(v) - GjCv) 

where all matrices have their elements In the quotient field R(v) formed from 
the ring R[v] of polynomials In the Indeterminate v with coefficients In the 
real ntmber field R. Our viewpoint on this equaton Is the free R[v]-module 
structure presented by Sain In 1975. Previous work has made clear how to answer 
questions such as whether G(v) can be proper and minimal. Herein we present 
some tools and partial results oriented to the question of whether G(v) can 
be stable. 

In the MacFarlane methodologies case, we have considered the problem of 
compensating a two-input, two-output, five-state Jet engine model In the class 
described by Michael and Farrar. Inputs chosen were main burner fuel flow and 
Jet exhaust area; outputs were thrust and high turbine Inlet temperature. 
Software was developed for constructing characteristic loci and a number of 
related parameters, and an effort was made to use the Ideas cf noninteraction, 
integrity, and accuracy to achieve an Introductory Jet engine control design. 
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The results are positive In naturei althou{;h we feel that more attention Is 
needed In this theory toward the actual design of reasonable compensators. 

Remark 

This work was essentially complete by February 29, 1976. The May 6 
date merely reflects delays induced by allowing time for thesis readers to 
provide rem£.rks and for scheduling of the thesis defense. 

Remark 

The main part of this work concerns the characteristic locus; the 
other result on polynomial methods relates to the previous section and con- 
umed a smaller time of the research. 


Remark 

Two portions of the work contained in this report are being prepared for 
presentation to the technical community. The first part, about assignment of 
the poles of the solution, is being prepared for the Fourteenth Allerton Con- 
ference on Circuit and System Theory. The second part, dealing with the 
characteristic loci studies, is being prepared for the Nineteenth ,.*dwest 
Symposium on Circuits and Systems. A very interesting feature of this pre- 
sentation, in the opinion of the authors, is the conjecture tljat difficulties 
in applying the basic method to jet engine control are intuitively due to 
interaction among the Inputs and outputs. There seems to be a possibility 
that such interactions can be related to the fact that certain jet engine 
characteristic loci do not separate into self-contained sub-loci as some- 
times occurs In less complicated examples. If this should indeed later 
turn out to be the cane, it would represent a substanti.il new Involvement 


K>,|| ,, 
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Cl topological concepts In feedback control systems, which have traditionally 
been more heavily Influenced by algebra. 

Remark 

The conclusions that more attention to compensator design Is needed in 
this theory, at least for the jet engine examples under study, arenot very 
prohibitive. New investigations to overcome this deficiency are now being 
planned. 


PART II 


GLOBAL NONLINEAR METHODS FOR JET ENGINE CONTROL 


n. Analog Simulation of Single Spool Turbojet 
T. C. Brennan and R. J. Leake 

b. Installation of DYNGEN Jet Engine Simulator 
J. C. Shearer and R. J. Leake 

c. Optimal Control Studies 
R. Basso and R. J. Leake 
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Ha. ANALOG SIMUI^MION OF SINGLE SPOOL TURBOJET 


This work represents an attempt to characterize the essential dynamical 
characteristics of a simple single-spool turbojet engine through simulation 
of low order system models on an Electronic Associates TR-48 analog computer. 
The objective was to gain insight into the most important dynamical constraints 
of such an engine, which can be applied to control studies of more advanced 
engines. 

Included is a derivation of a seventh order simplified simulation model, 
a derivation of an even simpler third order model, and simulation results 
from each. Details are given in reference (1) listed in the abstract of 
this report. 

The control problem studies is one of getting from "Windmill" (zero 
fuel flow equilibrium) to "Design" (a high thrust equilibrium) while taking 
into account surge margin and turbine inlet temperature constraints. Several 
control schemes were investigated. 

The principal constraints were imposed by the limited nonlinear equip- 
ment available on the TR-48 analog computer. The main simplifications involved 
the extensive use of linear approximations, single stage compressor dynamics, 
a linear compressor map, assumptions of a choked nozzle condition, and certain 
empirical relations based on design point equilibrium data available. Addi- 
tional simplification was achieved by limiting the model to a condition of 
20,000 ft. altitude at Mach number .8. We shall term this model the "Drone" 
system. 

The Drone System 

2 

Througljout the development, P stands for pressure in newtons/meter , 
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T stands for temperature in degrees Kelvin, («> represents mass flow rate In 

3 

kilograms/second, and p Is density In kllograms/meter . N Is rotor speed 
in revolutions /minute. Using a single stage compressor model the main dyna- 
mical equations are 


Auxiliary Relations 


Pj* - Rpjj- 

fA - “^B^A 
"a ■ S 


"3' “a 


Non-Llaear Functions 


Turbine 


"5 ■ 

Tj - fj(T^,P^,Pj,N) 


Nozzle 






Compressor 


Tj - f^(Tj'.N.P^) 




Normal l>;cd Seve n th Order Modol 

Using a number of assumptions detailed in reference (1) of the abstract 
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of this report, and normalizing about a selected design point the seventh 
order numerical model becomes: 


State Equations: 


A A 


dPj 

^ - 51.52 (^2 


dm, 

^ » 2188.27 (1 - PJ) 


^“3 » 8191.36 [1.05787P, - P. - .05787 

dt h 


^ » .93586 4];^ ^ 31.486i. + 21.435l,T, - 53.86 K 

dt £ 3 3^ 


dp— • A * 

^ » 37. 78^2 - 38.448 + .66849(L^ 


CP» Pj A A A 

dT ‘ *1-’’ r - VI 

D 


dt 


^ [3.12P, - 2.736lPeQl 

Pa 4 5 


.305N 


+ [.688013t:)2 - 1.0715 


Other Relations: 


“3 ■ “4 ^ ■ ‘>8^4 ’’2 ■ ^ 


Ta - 


» 



-18 



A 



** - ?5® 




T 3 - .64212 + .35788N 


Pj - 4.394N - 3. 394^2 
F - 0 {1.5486Pj - .5486} 


V- “ V. « V. “ A, « A. » .1 I ■ .0305 
14 5 14 


Normalized Third Order Simulation Model 
State equations: 


dP 

^ (.93586P^/Pg + 31.486) + 21.4353^1^ - 53.86 P^/p^ 


dp t« “ “ * 

•r-^ » 37.78u, - 38.448P, + .66849u^ 
dt 3 4 £ 


dN 1.258 2,. : „2, 

3F ~T~ ‘ ■ “3 ’ 

N 


Other Relations: 

All relations equations in the previous seventh order model arc valid, 


plus 
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A* ■ 

Pj-l 


W» • U). 


P30 


w, • 1.3009N - . 13982 (P, - 2 ^ 1777^. 

3 ^ . ^4 + .41688N - .0899P^N) 


Equilibrium Conditions 

Equilibrium conditions for the normalized seventh . order and third order 
models are the same. The ’’Design" equilibrium occurs vhen all normalized 
State variables are unity, and this corresponds to a design point of Mach .8 
and 20000 ft. Other equilibrium points must be calculated using a successive 
approximation procedure to solve the nonlinear equations which result by 
setting all derivatives zero. In particular, we define ’’Windmill" as the 

s ^ 

equilibrium which occurs when fuel flow *» 0 and nozzle area 6 ® 1 . 

The algorithm used to calculate equilibrium conditions is as follows: 


1. Set values of and 0 (the controls). 

2. As initial estimates, set 

- 1. N « = .5(ij + 1) 

3. Set in order 




'B 


1.3009N - ^ .«h88n2 - .0899^) 


.64212 + .3578SN^ 

4.394N - 3. 394^3 
(wj + 56.515 j3)/57. 515 

(53.86P^^ - .93586P^ij)/(21, 435^313 + 31.486ij) 
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*4 ■ 


r, - P*/e 


4. Set 


itk A A * 

r, /p u- 


and return lo 3, until convergence is achieved. 

This algorithm gave five place convergence in twenty iterations and showed 
a very nearly linear operating line relation of approximately 

Pj S 1.026362 - .0263 

A * * 

or equivalently (using the linear compressor relation “ 4.39AN - 3 . 39 A 62 ) 


«2 - .99405N - .0095 


Compressor Map and Surge Lines 

The equation used to determine the approximate compressor map of Figure 2 


is 

Pj/Pj - 5N/10^ - 2.94215023^2/^2 

where 

Oj- T2/288.3 *= 280.6/288.3 ® .97329 
«2 - P 2 /IO.I 325 X 10^ 


Using these relations and the equilibrium values 

P3^ - 28.076 X 10^ 

Pjg - 7.09 X 10^ 



I 
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Sb ■ 

- 3.48 X 10* 

ve obtain 

Pj - A.394P’N - 3.394i^ 

A 

This relation is graphed in Figure 2 for “ 1 along with the operating line 
of equilibrium points. 

We have chosen to define the surge line as a line running parallel to 
and above the operating line. This is rather arbitrary, but preserves the 
qualitative nature of the problem of getting from Windmill to Design without 
transgressing on the surge characteristic. Indicated on Figure 2 is a surge 
line specified by the relation 

Pj » 1 -j263w^ + .24105 
L inear Normalized Systems 

In order to estimate the dynamical modes of the various models, and to 
provide models for linear control studies, linear representations of the form 

** A6x + B6u 
at 

were determined, where A and 15 = -1^. Using the seventh order normalised 

equations and linearizing about the design equilibrium 

Xg - (1, 1, 1. 1, 1, 1, 1) 


wc have 
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* 

A 

n 

A 

«2 

*3 

• 

% 

A 

'2 

0 

51.52 

-51.52 

0 

0 

* 






‘2 

* ^ 

-2188.3 

0 

0 

0 

0 

> 


40159 

0 

-31967 

-8191.4 

0 

A 






% 

0 

0 

21.436 

-106.72 

52.893 

A 






A 

0 

0 

. 37.78 

-38.448 

0 


0 

0 

0 

61.97 

0 

A 

1 





N 

0 

2.2558 

-3.5131 

11.488 

-i:2587 


0 

0 

0 

0 

0 

61.97 

-8.9708 


A 

N 

0 

0 

40159 

15.343 

0 

0 

-2.5146 


B - 



0 

0 

0 

32.4219 

.66849 

0 

0 


e 

0 

0 

0 

0 

0 

61.97 

8.9708 
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The corresponding chlrd order normalized linearizations are 




-112.27 


-48.12 


A 



32.924 

0 


A 

N 

42.26 

47.444 


H 


2.8377 


-1.258 -4.096 


'B 


li 


•£ 

32.4219 

.66849 


These matrices were obtained tor the > zlues of V ■ A ■ and 1 
given on page 18. Other gain systers were studied by changing V, A, and I. 
A summary of digital computer calculations of the eigenvalues and dominant 
time constants for various gain systens is presented below. 


Gain System A . (V ■ .1, A ■ .1, I ■ .00305.) To obtain A and B matrices 
multiply bottom rows of all matrices by 10. 


7th order eigenvalues: 
(t ^ .0425) 


-31852. -141.43. -60.278 
-44.907 t j 329. 24 
-8.5923 ± j21.913 


3rd order eigenvalues: -89.288. -45.638. -18.304 

(t ^ .0546) 
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Cain Svfitctn B . (V - .1, A • .1, 1 ■ .0305) A and B matrices exactly as 

shown in equations. 

7th order eigenvalues: -31892, -88.732, -61.218 

(t ^ .296) *25.20A, -3.3818 

-33.738 ± J333.68 

3rd order eigenvalues: -81.219, -32.309, -2.8386 

(t J .352) . . 

Gain System C . (V - .01, A - .1, I ■ .0305) to obtain A and B matrices, 
multiply all rsws except the bottom one by 10. 


7th order eigenvalues: -318962, -81A.82, -618. 5A 

(t .339) -310.44, -2.9529 

-325.38 ± j3339.7 


3rd order eigenvalues: -807.68, -316.2, -2.9166 

(t y .343) 


Gain System D . (V - .01, A ■ .1, I ■ .0305) to obtain A and B matrices, 
multiply all upper rows by 10 and the bottom one by 1/10. 


7th order eigenvalues: 
(t y 3.423) 


3rd order eigenvalues: 
(t y 3.420) 


-318966, -806.7, -619.58 
-316.34, -.29218 
-324.17 ± j3340 

-807.25, -315.57, -.2924 


Simulation 

The actual simulation was conducted at the University of Notre Dame's 
Analog Computer Laboratory. An Electronic Associates TR-48, and n few small 
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TR-20 and TR-10 analog comptucrs were used. Patching diagrams were obtained 
with the aid of ANSIR 3, a digital computer program which generates all perti- 
nent analog simulation information, given the differential equations. 

Control Problem 

As stated earlier, the control problem is to schedule fuel In order to 
accelerate the engine operating state from Windmill to Design equilibrium con- 
ditiono. However, maximum acceleration potential is limited by compressor 
surge or stall. As fuel is increased to accelerate the engine, pressure ratio, 

A 

P^/P^, increases and airflow, decreases. When stall is reached the pres- 
sure will drop and will cause a similar decrease in compressor efficiency. 

It is desirable to design controls which will cause die compressor to approach 
the stall line, however, since these controls tend to give a faster engine 
response. 

Four different fuel control concepts were investigated on the simulation. 
They consist of an open loop control, !■* -lear feedback design, and two nonlinear 
controls developed by empirical methods. Their ffects on engine performance 
and surge constraint were indicated by time response curves and trajectory 

plots. 

Control 1 . Here we use the term open loop to refer to this control. It 
simply involves letting 

* 

Wj = 1 

A 

This corresponds to the value of at design equilibrium. 

Control 2 . Control 2 fuel scheduling is a linear combination of compressor 
and turbine pressures given by 
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Sj - .99329 Pj + .00608 


Control 3 . This control Is the product of compressor airflow and rotor 
speed. It Is defined as * 

* • c§ 

"f ■ “3 “ 


Control 4 . This nonlinear control Is given by 


^ A O A 


These controls were studied for both the seventh order and third order 
models given above and detailed results are given In reference (1) listed 
in the abstract. 

The exercise was generally successful, yielding valuable experience in 
determining the essential dynamical characteristics of jet engines. The Sj.rau- 
lation results were reasonably representative of the types of response one 
would obtain using more sophisticated models. 


lib. INSTALUTION OF DYNGEN 


This study Is concerned with the setup and utilization of DYNGEN, a 
FORTRAN program for analyzing steady state and transient performance of jet 
engines. DYNGEN is based on three earlier computer programs; SMOTE, GENENG, 
and GENENG II. These programs allowed the user to analyze the steady state 
behavior of severa] different engine configurations. However, as the need 
to predict the transient performance of these engines became more Important 
in preliminary design and control studies, DYNGEN was developed to provide 
this added capability. The result is a digital program which can perform both 
steady state and transient calculations on a wide variety of engine configura- 
tions without any significant modifications to the program Itself. 

Because of DYNGEN ’s size (slightly less than 32,000 words of computer 
storage) certain decisions had to be made concerning the storage, editing and 
running the program. It was expected that DYNGEN would be used quite often 
in subsequent control studies. This made the reading and compiling of its 
5,000 plus cards too cumbersome and costly. Furthermore, the control sub- 
routines FCNTRL, N02CTR, and DISTRB would have to be edited, as new control 
schemes were devised. Again, rereading and recompiling the entire program 
for each new control would prove nonfeaslble. To overcome these and other 
difficulties, the Notre Dame PURT (Program Utility Routines for Tape) system 
available on the IBM SYSTEM/370 MODEL 158 was utilized. A description of user 
options and the JCL (Job Control Language) required to initiate them can be 
found In more detail in reference (2) listed in the abstract of this report. 

Once DYNGEN was made readily accessible, a manual providing the Informa- 
tion required to run the program was written. Because It was assumed that only 
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one engine configuration would be used In further studies, the Information 
In the original DYNGEN report cov’ld be vastly reduced and simplified. As 
a result, the time and effort required of a beginner to learn the particulars 
of the program were kept to a minimum. 

Tape Initialization 

The Initial decision as to how to store and access DYNGEN was rather 
crucial. Ease of handling and editing and minimizing compiling and execution 
times were Important considerations. Cards proved too slow, cumbersome, and 
costly to be used effectively. Disk storage is fast and convenient, but the 
storage requirements of the program were rather excessive. Magnetic tape 
seemed to be the best alternative. Tapes are reasonably fast, cheap to use 
and, with the available PURT system, quite easily accessed and edited. 

The DYNGEN Program was received from NASA on a seven- track tape with a 
556 bpi density. This had to be transferred onto a more commonly used nine- 
track, 1600 bpi tape. After securing a 600 foot tape from the Computing Center 
(tape number M0151, code name RAT) and submitting the NASA seven-track tape to 
the operator, the following JCL program transfers the information from the seven- 
track tape to the new nine- track tape. 


//1AI>F JOU <CF .ft) «STUtlCUTID.TIMt = 09,MSGLEVtL=a.l) 

/♦HDU1L t'Kl II Luf AL 
/t JOI'I’AUM 

/•M llil' iniloLIJ.UlJiPTM 

// lArc pon=irc>c.i.iiLK 
//.SY;.t'KPil Uii i>1.v)ul = A 
//StSiili I»U = . 111)17 =1 AI'L7 . 

// l)L*»: ll M.: 1 . ll' li icL I , klCI M , UJl Ct = 00 . OLKS 1 ZC = 00 0 ) » V0L = 5C K = 0t NCtN « 
// L)ii'iL = ( 1 .itLt 1 .i)i;.r={(ii Dii’Asn 1 
//S7M'7/ UU W.Kl |.t>) , 

// l.tifiLL - I 1 . NL I . vl»L = M K = nu 1 M , UNI 7 -1 Al'L » 

// uC LF ■•’.il li.Uil.CL:t;l).oLKblZL-aOU) 

//St Ml) O'J • 
i'.UJl.«A7t NAxFLO.'.rl 
htCOMU FlLLUsUWi l,Ht ) 

/« 

// 


Of & fc- 
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Once the new nine- track tape was prepared, the entire FORTRAN program 
was punched out on cards. This was to make correcting easier and to provide 
a backup copy in the event the tapes were accidently misplaced or destroyed. 
The following JCL program performs the punching operation. 

//PlIf.r.H JOb ( CF t L0b9> .STUOLNTiDf T1KL = 01 .HSl.LL^/ELr I J, 1 1 

/♦JOUt'ARf'i tAt<US = bUO0 

/•ROUU PUNCH LUC AC 

/♦SHOP HuJbl.KAT 

//SltMl E*f.C POS=It.i>PTPCH 

//StSMKiUT Ut) bTS()UT = A 

//SrSUTl l)|| USH = uYr.UUJ,UlSP=(OLO.PASS) ,LiUCL=(l.SL».VOL = SCR=M0151, 

// UHITsTaPl 
//STSUT2 DO srsOUTsb 
//STSIU t)b * 

POUCH TYPOKOsHS 
/• 

// 

The resulting source deck was divided into four separate source modules. 
The first three modules consist of the three control subroutines: FCNTRL, 

NOZCTR and DISTRB. The remaining source module contains the rest of the 
subroutines, dummy mainline, and control functions. Each of these source 
modules is read in turn into the compiler and translated into four sets of 
machine Instructions called object modules. In order for these object modules 
to be executed, they must be processed by the loader. The loader accepts all 
the object modules from core, resolves the external references, searches the 
subroutine library, and combines everything into an executable form. Con- 
trol then passes to the program for execution. 

The advantage of segmenting DYNGEN in this manner is straightforward. 

If any of the source modules has to be corrected or replaced, only that 
module has to be recompiled. The new object code from the edited source 
module can then be combined with the remaining old object modules by the 
loader and executed. Because the bulk of DYNGEN is put into one module 
which should never need revision, the savings in time and computer costs is 
quite significant. Another added feature of this system is that new source 



31 


and object modules can be loaded and executed without destroying the old 
nodules. The programmer then has the option of which modules to save. 

The following is the JCL program needed to segment the entire DYNCEN 
program into the four source modules, compile each, and store both the 
source and object modules back on the tape. It should be noted that this 
program should only be used to segment the original DYNGEN. 


//OHL'AKUP JUo .KotJl) .i.TUOLwnu.TU'.LsU^J.NSbl-LVtLsil.ll .KE&I01J=19;?K 
/• tiWtJsSiO 

/•Stflil* Nuljl.KAT 

//sa»'i extt 

//SUHa CXLL »‘UKU.J.T;U't-M0l3l 
//STU'A EaLC l'UKIi,uUf« 

//SUuKCE.SIbJiM Ely ♦ 

• / AJU UAHL = iJTKbLUl 

./ IULK=1U0 

» » 

• * tr.ltU *.Lu t-imus tXCLPT Tut » 

* IHKLL LOi-jIKUL EUbKUUUWLS. * 

♦ * 


//STLl‘1 E>.lC KiU 1 Kol<li.rgA,lL=UYtiGEEl 
//STLP'j txtt Pul'liswUK 
//S 0 Ut<CL. 3 ti,l;'j Uo ♦ 

,/ Atu NAilt = iiJ Mull 

• / l.u^iuEK la. t 

♦ ' • 

« IIIL UlSlUU SUbKUUT IM. . * 

« ♦ 

//STtPG txii; h'LiMI-ul<H.i,A.lL = UlS1i{b 
//S1U’7 txLC l''->hlbyjbl4 

• / AUU UAnL = l till Kl 

• / l 4 U-‘-ULl< 14L A 1 = Jl u U I 1 r.bH- X UO 

4 . 44 «* 4 »«** 44 »»* 4444 i 4 >« 444 . 4»»»4 

• LhTlU IIiL FlMKL SUbuOUllUL 


//STt.»T; rX5-t I'Oi, lKoliM,i.Arit. = hCNll<L 
//SYLi'? r.xLt CbN ( Sbljll 
,/ AUil MU41C 1 1< 

,/ IgUf'ULii UL A 1 = 1 ly 0 I I 1 ill) 


LtgiLK IHL r.u^ClI. bOiiKCUl lUtl. 


4 %. 

••• 


//STU'lO LjU. 1. t'UMHiKii,r;i.^L = i.’U/CH: 
//SlLMl LxLl I'Uulubl . 1 

// 


The tape is now in its final form. It can be executed, edited or both 
with just a few simple JCL cards. Each option the user desires rciiuircs one 
or more JCL curds in "block" form. Fnrtticrmorc, each desired option may 
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require several preliminary JCL "blocks" for Inltlaliicatlon - a tape cannot 
be edited until It has first been read Into the main core memory. It Is 
Important, then, for the user to understand how each "block" works as well 
as be able to reproduce on cards the JCL programs that follow. 

Program Initialization 

Before anything can be done with the tape, it must first be put on the 
tape drives and read into core. This Is accomplished by the following JCL 
cards . 


//OYNGEN JOB (CF .F081 ) ,STUDE;NT10iTlMEs02«MSGl.tVEL= ( OtO ) |RE1G10N=320K 

/♦JOBPARM L1U'ES=20 

/•SETUP K0151.RAT 

//STEPl CXEC PURTINIT 

//STLP2 EXLC PURTl.N.TAPt=M0151 


The emphasis here will be on how to execute DYNGEN and how to edit the control 
subroutines. If the main body of subroutines, located in the module named 
DYNGENl is to be revised, TIME on the JOB card should be set to at least 
three minutes for editing purposes alone. 

Editing 

Once the source and object modules are loaded into core, they can be 
easily edited. \'Jhcn DYNGEN was divided into four parts, each module had its 
statements numbered in Increments at 100 starting with 100. By specifying 
any individual statement number, any statement can be changed or deleted with- 
out affecting the rest of the source module. Also, the entire source module 
can be changed by submitting a new dock of source statements. Regardless 
of the editing method chosen, each revised module must be recompiled and. If 



it is to be saved, read back onto the tape. To get a text of all statements 
and their assigned sequence numbers, refer to the latest compiled listing or 
have the necessary modules printed out. 

There are two print options available to obtain a listing. First, a 
listing of everything on the tape can be obtained by the card: 

' //SUPM.EXtC PUKTLIST 

I 

If only Individual modules listings are desired, the following JCL cards 
are required. 


//S1LP4 ExLC PUMLIST 
//PUK^tlST .SToirj LU * 

« 

• LMtH THL KAHL tr TUl LOAD Mlit/Utf.S 
» To HL USILl) UH COM»ttUnVL tAivtSi 

• SlftnliKti lU IhL HKi 1 tOLU.SIj 




Once the sequence numbers of the statements to be edited are known, it 
is a simple matter to make revisions. First, the corrected FORTRAN statements 
should be punched on cards with the appropriate sequence numbers in columns 
73 to 80. If a statement is to be inserted betweer'. two other statements, the 
sequence number on the new card must be some number between the sequence 
numbers of the other two statements. The following JCL cards are required 
to make statement changes. 

ixcc puKisuurt 

//SOUKCC .:.YS1N to • 

,/ CHAtJC-t UA.lLslLinLK NOUULt t/.NC » 


tt.lLN ALL M.W STAUflLtjTS 
COLUMNS • 


4 
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The following JCL cards are needed to delete a statement or series of 


statements from a module. 


//STtl'7 Eiktc PUKlSuUU 
//SOUUCC.STblN i;u • 

• / ULLLIt i/Lwl-lvi&wtcsn 

/• 


This will delete statements n through m. If only one statement Is 


to be deleted, let n - m. 


Sometimes, it is more coc.venitut to resubmit the entire revised module 


rather than edit each individual statement. To change an entire module, a 


new source deck is punched out on cards and is entered in the following 


manner: 


//STtPS ZXtlC l»Ur<TSOUR 
//SOUUCE.SYi.Ii'. DU • 

«/ RLPL I c-.'.- I lK i'iOOULl I'iAi'E) 

•/ UUKUCK Nlw 1=1U0<1 NlK= 1C0 

» enter SOUHCE ROUUI.L UtCK. 

* 

/• 


To compile a newly edited source module, the following statement is required. 


//STEPO EXLC PUf<lFURH,NA.'iE=<EraEK P.OUULt fJAKL I 


If more than one modulo has been revised, the above card must be supplied for 


each. If the DYNGENl module is being edited, an extra card may be necessary. 


This is to insure that there is sufficient disk space available for the com- 


pilation. 


//SUPt fXLC PUHTFU«M.JjAr.r=DYI.&Lfjl 
//» um .STUU•^^ uu SKACLMtlL.UtiJ) 
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To execute DYNGEN after all source and object modules arc in their 
desired form requires the following cards: 

//STCP9 LXte PUHTFG0.(;At'.c;s0tfJ5ENl 

//CO.KTOof UUl OD OUPsii^c-imHAsS) , SPACLs I TRK « (10 • 1 0 > ) ilMT = STi.DA 
//CO.&TSIw UJ • 

• • 4 » 4 **«»**»«***«******'********»*** 

• CNICK DYNGEN INPUT CAROS. • 

• « 

/• 

Control then goes to the dummy mainline program in the DYNGENl module. 

The editing techniques explained thus far change only the modules as 
they appear in core. To transfer the corrections back to the tape is the 
function of the PURTOUT card below. 

//STEP 10 EXCC PUKTOUT.TAPESK0151 

It is important to remember that if this card is not used, what was origi- 
nally on the tape will remain, even if changes had been made to the source 
and object modules. Also, when more than one option is used during a parti- 
cular run, the /* should be omitted between blocks. It should, however, be 
used after the last block. 

Further discussion on the Installation of DYNGEN, together with some 
simple control studies is contained in reference (2) listed in the abstract 
of this report. 


lie. OPTIMAL CONTROL STUDIES 


The analog and digital simulation sutdies of the previous two sections 
are aimed at establishing models of jet engines for the purpose of experi- 
mentation with various control system designs. We are concerned with both 
linear multivariable and global nonlinear approaches as alternatives to 
linear quadratic Rlccatl methods. Here, we describe some preliminary studies 

by R. Basso and R. J. Leake on the global nonlinear approach. The work is 

« 

incomplete as of the end of the grant period and will be completed in the 
summer of 1976. Some details are given In reference (6) listed in the abstract 
of this report, and will be presented at the 1976 Joint Automatic Control 
Conferences. 

The linear multivariable approach 's of the previous sections a-*m at 
various criteria common to good regulation about design equilibrium points. 

They depend on scheduled nominal inputs for large excursions which are often 
part of the engine design itself. A global nonlinear approach actually takes 
scheduling and local control action into account simultaneously by specifying 
the total control, and in this manner, overlaps the control and engine design 
problems. 

As a preliminary global study, the problem is taken as one of going from 
idle to a high thrust design equilibrium point in minimum time, while observing 
side constraints on surge margin and high turbine inlet temperature. 

The usual approaches of Optimal Control have been considered, but only 
within the more modern, more systematic, and more general framework of Mathe- 
matical Programming. If side constraints are included through penalty and 
barrier methods, the control problem can be viewed as an uncons tr ained mini- 
mization of an objective function J(u) where 
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u(0) 

u(l) 


and 


N-1 

J(u) •= K(x(N)) + I L(x(t),u(t),t) 
t-0 


with the differential equation approximation 


x(t+l) - x(t) + f (x(t),u(t)). 

Standard algorithms such as the Fletcher-Reevcs Conjugate Gradient and 
the iix»dlfled Davldson-Flctcher-Povell quasi-Newton algorithm with self-scaling 
and automatic restart can be applied directly to this problem with J(^) being 
evaluated for t ■ 0, 1,..., N-1. The gradient 


u ‘ u(0) u(l) u(N-l) •' 


is determined by solving the adjoint equations 


X(t) - X(t+1) + X(t+1) V^f(x(t), u(t) + V^L(x(t), u(t).t) 

X(N) » V^K(x(N)) 

f 

for t • N - 1,..., 2, 1 and observing that 

Vy(t)J - ^(t+l) V^axCt), u(t)) + V^L(N(t), u(t),t) 

Each of the algorithms mentioned above requires a line search for minimization 


h(a) - J(u+ad), h'(0) - V^J(u) 
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along directions (usually not the gradient direction) determined by 
the particular algorithm. A cubic approximation of h(a) depending on h(0), 
h*(0) and two other values, h(a), h(b), as suggested by Lasdon, seems to be 
an efficient basis for carrying out a fairly simple line search in one step. 
Refinements of this approach are being sought, as it is essential that the 
number of evaluations of J(u) and he kept to a minimum due to the 

amount of computation required. , - * 

If the variables are normalized about a hlfh thrust design point and cer- 
tain simplification^ are assumed for the nonlinear functions involved, a com- 
plete numerical model which has the essential characteristics of the engine 
is given in Section Ila of this report. 

The single control variable is WFB and the problem is to drive the system 
from idle (actually windmill here with WFB “ 0) 

pj - 1.774 - .5384 N^, - .5461 


to design (with WFB ■ 1) and 


A 



1 


A 


P - 1 
4 



1 


in minimum time subject to surge margin and turbine inlet temperature con- 
straints approximated respectively by 


P. < 1.25 
4 — C 

P^ < 1.25 Pg. 


Notice that a good second order approximation to the problem can be 
obtained by taking as a control variable. Once the optimal is obtained, 
It Is used to determine the WFB necessary to yield P^. This technique was 
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used to obtain a Dynamic Programming solution to the problem with time step 
oitc 7t ■ .002. Generally, there arc three regions of control: 

(1) rides the surge margin constraint boundary (P^ • 1.25 N^). 

(2) switches and rides the inlet temperature constraint boundary 
(P^ - 1.25 Pg). 

(3) P^ leaves the inlet temperature constraint boundary and adjust to 
get all states to final design equilibrium value of unity. 

The optimal time irom idle to design is slightly less than one second. 

A preliminary open loop solution, using a conjugate gradient algorithm 
with time variable penalty functions to account for state constraints and 
terminal condition, has been obtained but there are some constraint viola- 
tions and the solutions are, as yet, unsatisfactory. 

A 

The constraint violations which occur in P, are due to exterior penalty 
functions. Feasible c ction and barrier approaches are being studied to 
remedy this effect. 

The fact that the problem Is time optimal and the optimal number of steps 
N is initially unknown docs not appear to offer any real difficulty as a few 
rough computational runs quickly give one a good estimate of N. 

Future studies will involve experimentation with various conjugate gradi- 
ent and quasi-Newton algoritluns in combination with penalty, barrier, and aug- 
mented Lagrangian techniques. Also, it appears that the state constraints can 
be taken to be of very simple form (such as Xj^(t) ^ 0) and this should open 
the way for the development of simple feasible direction methods to handle 


the hard side constraints. 
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ONGOING AND RELATED ACTIVITIES 


In addition to the activities described in the preceding two sections, 


other grant activities, ongoing and related projects, are taking place. 


1976 Joint Automatic Control Conference 


With the assistance of J. Zeller and E. Lehtinen of the NASA Lewis 
Research Center, the investigators have arranged a special sponsored session 
at the Joint Automatic Control Conference, Purdue University, Lafayette, 
Indiana, July 28, 1976. A brief description of the session follows: 

Application Session 2 


Theme — Application of Control 

Title "Jet Engine Control" 

Chairman M. K. Sain, University of Notre Dame 


Vice-Chairman — B. Lehtinen, NASA i.ewish Research Center 
Paper No. 1 

"The Role of Multivariable Control Techniques in the Design of 
Turbine Engine Control Systems" 

C. Skira, Aero Propulsion Laboratory 
Wright Patterson Air Force Base, Ohio 

Paper No . 2 

"Multivariable Control Design Principles with Application to the FIDO 
Turbo fan Engine" 

R. L. DeHoff and U. E. Hall, Systems Control, Inc. 

Palo Alto, California. 

Paper No. 3 

"Simulation of a Turbofan Engine for Evaluation of Multivariable 
Optimal Control Concepts" 

K. Seldner, NASA Lewish Research Center 
Cleveland, Ohio 

Paper No. 4 

"Application of Multivariable Optimal Control Techniques to a 
Variable Area Turbine Engine" 

E. C. Beattie and W. R, Sprock, Pratt & Wliitney Aircraft, Inc. 

East Hartford, Connecticut 
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Paper No. 5 

"AlternaMvc Methods for the Design of Jet Engine Control Systems” 

M. K. Saint R. J. Leake, R, Basso, R. Gejji, A. Maloney, and 
V. Seshndrl, University of Notre Dame 

Notre Dame, Indiana 

This session is officially sponsored by the Control Systems Society 
of the Institute of Electrical and Electronics Engineers. 

Dominance and Inverse Nyquist Studies 

Mr. A. Maloney, a teaching assistant not receiving direct support 
under this grant, has been developing software to help in the evaluation of 
the methods of H. Rosenbrock as applied to jet engine models. Selected 
curves from this work have been used in Paper No. 5 of the JACC session 
described above. The preliminary work is expected to be complete during 
the coming sununer. • ■ ■ 

Classroom Evaluation 

Jet engine models were used extensively as project models in the course 
EE 555, Multivariable Control Systems, at the University of Notre Dame in 
Spring 1976. Several insights into various frequency domain techniques were 
obtained in this way, and one or two further projects may possibly grow out 
of this class work. 

In the course EG 551, Mathematical Programming, various algorithms were 
Investigated for computation of time optimal controls for the simple jot 
engine model discussed in Parts lla and lie of this report. 


